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ABSTRACT 

EGRET (20 MeV to 30 GcV) on board the Compton GRO has observed 
high-energy emission from about ~ 40 — 50 Active Galactic Nuclei. Theoretical 
models of this emission based on the upscattering of thermal disk photons 
by coohng, relativistic electrons can successfully account for the EGRET 
observations, but they predict a considerably greater X-ray flux than that 
actually observed in a majority of these sources. This inconsistency may be 
an indication that the particles are energized during the Compton scattering 
process, since the X-ray emission is produced by the lowest energy electrons, 
whose density may be relatively small due to the acceleration. Such a situation 
may arise as a result of resistive field generation in electromagnetic acceleration 
schemes, which we here explore. A key feature of this model is the assumed 
existence of a current associated with the azimuthal component of the 
underlying magnetic field by a slight imbalance in the energy distributions 
of outwardly moving, relativistic electrons and protons produced at the disk 
surface via shock acceleration. The generation of an electric field (via magnetic 
field line reconnection) is thus required to maintain the current in the presence 
of a resistivity induced by the radiative drag on the relativistic electrons. We 
show that the resulting spectrum can exhibit a significant deficit of X-rays 
compared with 7-rays. In addition, due to the uni-directional flow of the current 
associated with S^, this model would predict that the electrons are energized 
relative to the protons in the outflow only on one side of the disk. They should 
be decelerated on the reverse side. As such, we would anticipate that any given 
blazar should have a ~ 50% probability of being 7-bright, which appears to be 
consistent with the observed ratio. 
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1. Introduction 

Many active galactic nuclei (AGNs) produce radio jets or strong radio core emission 
whose origin is still not completely understood. Because many quasars are also luminous 
X-ray sources whose strength is correlated with the radio core emission (Browne & Murphy 
1987), it is believed that understanding the high-energy emission, which probably originates 
close to the central, supermassive black hole, will yield important insights into the formation 
of these jets and the radio emission mechanism. During the past few years, EGRET on 
the Gompton GRO (designed to measure photon energies from 20 MeV to 30 GeV) has 
detected about 50 AGNs, of which ~ 40 are high confidence detections (Mukherjee et 
al. 1997). Most of them are quasars, at least six are BL Lac objects, and at least 13 
are Optically Violent Variables. The vast majority are radio loud, flat spectrum objects 
ranging from a redshift z of 0.03 to more than 2.0. Typical blazar (i.e., AGNs whose 
dominant observational characteristics appear to be due to beaming effects) spectra suggest 
the presence of a break below ~ 10 MeV if the 7-ray component is to connect smoothly 
with the observed X-ray spectrum at lower energies. This gamma-ray emission is most 
likely associated with a particle acceleration region near the black hole, thereby providing 
tantalizing clues into the physical processes acting in this environment. 

A wide variety of particle energizing mechanisms have been proposed for the jet 
formation, ranging from radiation pressure-driven flows (Abramowicz & Piran 1980) to 
electromagnetically-driven self-similar jets (Li, Ghiueh & Begelman 1992). The most 
enduring models are the electromagnetic ones because they alone appear to have the 
capability of accelerating particles to the relativistic speeds necessary to account for the 
radiative characteristics of superluminal sources. The retardation of the particles due to 
the photon Gomptonization in the intense radiation fleld of these compact regions has 
been recognized ever since it was discovered that high-velocity jets could be decelerated by 
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particle-photon collisions. This mechanism is sometimes used to account for the low terminal 
Lorentz factors (700 ^ 10) observed in superluminal sources (Phinney 1982). Subsequently, 
Melia & Konigl (1989) demonstrated that the Compton deceleration due to scatterings with 
accretion disk photons could transfer most of the jet's mechanical energy into an observable 
flux of high-energy radiation, accounting for both the low values of 700 and the relatively fiat 
X-ray spectra seen in some sources. Indeed, inverse Compton scattering models involving 
a cooling population of leptons have been successful in reproducing many of the EGRET 
observations. However, the incomplete coohng of these leptons in the presence of the disk's 
nonuniform photon field results in spectral breaks between the 7-ray and X-ray regimes of 
~ 0.5 (Dermer & Schlickeiser 1993), a result that does not account for the large variety of 
spectral breaks seen in the EGRET AGNs (von Montigny et al. 1995a). While the presence 
of additional energy-loss mechanisms can increase the spectral break to ~ 0.7 (Dermer, 
Sturner & Schlickeiser 1997), it seems that the underlying problem with standard cooling 
models is the relative overabundance of low-energy electrons in the assumed power-law 
distributions. To this end, recent theoretical efforts invoking Comptonization have assumed 
the presence of an electron acceleration process to shift all of the leptons to high energy, 
thereby producing a low-energy cutoff in the electron distribution. This energizing is 
assumed to occur either before injection into the emission region (e.g., Bottcher, Mause & 
Schlickeiser 1997; Dermer, Sturner & Schhckeiser 1997; Sikora et al. 1997) or throughout 
the interaction zone (Marcowith, Henri & Pelletier 1995; Bednarek, Kirk & Mastichiadis 
1996). While these models have had success in accounting for both the X-ray and 7-ray 
data, the physics of how these conditions arise remains unsettled, particularly with regard 
to the specifics of the acceleration scheme in the latter class of models. 

In this paper, we will be addressing two principal issues. First, we attempt to develop 
a physically self-consistent picture for the mechanism that produces the required particle 
distribution with a low-energy cutoff. Second, we will follow the evolution of the particle 
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properties with distance away from the source of injection, and we will determine the effect 
of this on the cumulative photon spectrum. The view studied here is that low-energy 
electrons are first energized in a thin transition region in order to produce the current 
associated with the presence of a curled magnetic field and they then continue to be 
accelerated during the Comptonization process in order to maintain this current. Unlike 
previous models of this type, the acceleration scheme is well-specified and the electron 
distribution is determined self-consistently. The overall process is assumed to occur within 
an approximately force-free magnetosphere, as first described by Blandford & Znajek 
(1977). Following Meha & Konigl (1989), we assume the presence of a luminous accretion 
disk (often invoked to account for the observed "blue" bump) whose radiation scatters 
with the accelerating particles. The magnetosphere thus retains a small resistivity due 
to the particle-photon interactions that couple the acceleration of the particles and their 
concurrent radiative emission, but the fields remain close to their force-free intensities. This 
resistivity induces a component of the electric field parallel to the underlying magnetic field 
that drives the outflow. The object of our model is to determine the individual particle 
dynamics and the upscattered radiation spectrum self-consistently. The merit of such an 
approach is evident, in that the energized photon distribution depends critically on the 
particle energies. In §2, we describe a scenario with which the magnetospheric current may 
be produced, and we discuss how this current would be a source of 7-ray emission. Our 
calculations and results are presented in §3, where we also discuss their relevance to the 
observations. 

2. The Pcirticle Dynamics and High-Energy Emission 

For simplicity, the disk is taken to be a geometrically flat, blackbody emitter with 
inner radius Rin, outer radius Rout and a radius-dependent surface temperature as given by 
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(Pringle 1981) 



T(R) — Tjnax \ 7 



49 Rr ^ ^ 



36 R 
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D \ 1/2 
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(1) 



where the maximum temperature T^^ax = {^GM^M /bdnaBRmaxY^'^ (where is the black 
hole mass, M is the rate of mass transfer through the disk and is the Stefan-Boltzmann 
constant) is attained at the (equatorial) radius Rmax — (49/36)i?i„. The system is naturally 
described by a cylindrical coordinate system with an origin centered at the black hole and 
with the z direction normal to the disk. In this first set of simulations, the magnetic field B 
threading the disk is assumed to be axisymmetric and, as a result of disk rotation, to have 
an azimuthal component S^. 

The gravitational energy dissipation power is characterized by Eaccr = GM^M / Ri^. 
About half of this energy is dissipated throughout the disk. In the picture developed here, 
we assume that energy is transferred to a sub-population of electrons and protons via shock 
acceleration at a rate rjEaccr throughout an active region that extends from the inner edge of 
the disk out to a radius Ra "C Rout {v characterizing the efficiency of the particle energizing 
process). Using the above expressions for T^ax and Rmax, the rate at which the particles 
gain mechanical energy is expressed as 

EmecH = 4.2 X 10^^ ergs s^ (^) ( ) ' ( Vi^V , (2) 

""^"^ ^ VO.I/ \5xWkJ \W MqJ \3RgJ 

where Rg — 2GM^/c^ is the Schwarzschild radius. 

As a result of shock acceleration, the energized protons are expected to acquire a 
relativistic energy distribution characterized by a power-law in the Lorentz factor 7. While 
the electrons are also expected to acquire a relativistic energy distribution, the details of 
this process are not well understood. As such, we make the reasonable assumption that 
the energized electrons are also characterized by a power-law distribution in 7 but we do 
not require that the spectral indices for the two populations be equal. In principle, a 
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complex magnetic field structure would impose a complicated angular distribution for the 
charges. We note, however, that the observation of GeV photons in the EGRET spectrum, 
coupled with high UV/X-ray luminosities and rapid variability, drastically constrain the 
beam "size" of the high-energy emission (e.g., von Montigny et al. 1995a). This in turn 
limits the range of possible magnetic field geometries (i.e., B^<^ B), and suggests that the 
dynamics of the outflowing particles may be approximated by a more or less unidirectional 
motion. (The resulting 7-ray emission does, however, depend on the angular distribution 
of the electrons. This point will be addressed below.) For simplicity, we further assume 
here that = and take the energized particles to move directly outward along the z 
direction (since i?<^ <C B). For a uniform active region, the initial relativistic electron 
(proton) distribution function is therefore ne(p){'^'iZ = 0, r) = C 7~°e(p) for 1 < 7 < 7„;e(p), 
where ne(p)(7; r) d'j is the number density of electrons (protons) at (r, z, (f)) with a Lorentz 
factor between 7 and 7 + d'-f. Note that since Trip >> rUe, almost all of the mechanical 
energy is carried by the protons for comparable values of spectral indices. For simplicity, we 
assume throughout this paper that ccp = 2.5. The results of our calculations will therefore 
be insensitive to the value of the upper bound for the proton distribution •ju-p- The value 
of the upper bound ■yu;e for the electron distribution is determined by balancing the rate of 
acceleration within the shock region with the combined rate of cooling due to synchrotron 
and Compton processes (Begelman, Rudak & Sikora 1990). It is important to note that 
synchrotron cooling dominates over Compton cooling within the shock region since the 
magnetic field is expected to be highly disordered within that region (thereby resulting in 
large pitch angles between the local magnetic field direction and the electron motion) and 
the magnetic energy density is appreciably greater than the photon energy density. This 
therefore leads to a maximum Lorentz factor 7^ = 1.16 x 10'' (B/10'^G)~^^^. For simplicity, 
we adopt the fiducial value 7„ = 10^ in our simulations. In contrast, since the magnetic 
field becomes ordered away from the disk and the electron motion becomes nearly parallel 
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to the field direction, Compton cooling processes dominate the dynamics of the fiow above 
the disk. For a uniform active region, the value of the number density of energized particles 
is found by equating (2) with the expression for the mechanical energy flux carried by the 
protons Ep (7p)mpC^27ri?^cneo, where (7^) is the average Lorentz factor for the protons, 
rrip is the proton rest mass, and the factor of 2 takes into account both sides of the disk. 
Assuming that the flow is charge neutral, the number density of electrons energized by the 
shocks can therefore be expressed as 

(We note that for the assumed power-law index ctp = 2.5 for the proton distribution, 
(7p) = 3 so long as 7„;p >> 1.) 

The shock mechanism which energizes the electrons and protons acts independently of 
the large scale magnetic field structure through which the outfiowing particles then move. 
We therefore assume that the particles must undergo an adjustment in their distribution 
within a (presumably narrow) transition layer in order to produce the current associated 
with the presence of i3<^. We do not attempt to develop a complete account of this 
process here but rather postulate that it occurs within a height Zq due to the action of an 
electric field generated by magnetic field line reconnection so as to attain a self-consistent 
magneto-dynamic system. This assumption is reasonable as long as the particle energy 
density u ~ Ernech/'^R'ic is much smaller than ub = B^/Sir, thereby constraining the 
radius of the active region to a value Ra » {%Emech/ B'^cY^'^ ■ A comparison between rieo in 
Equation (3) and the critical number density necessary to produce the required current 

, Z . ^ . 0.56 cm- (?±\ (4) 

' ec 47rei?<^ \WG) \Rg ) \WMq) ^ ' 

(where Rif, is the characteristic radius of curvature of B^), then suggests that ric is probably 
due to a slight imbalance between the electron and proton velocity distributions since 
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rtc <^ rieo- Throughout the transition zone, the electron distribution undergoes either an 
upward or downward shift, depending on the direction of B^. Since the electric field shifts 
all the electrons over a given distance by the same A7 (ignoring pile-up effects at 7 = 1 
for a downward shift; these electrons do not contribute to the 7-ray emission), the electron 
power-law distribution is maintained throughout the transition zone but the value of the 
lower cut-off 7^ may be raised. As such, we shall assume a power-law for ne(7; z^^r) at the 
transition zone boundary and take 7/ to be a free parameter characterizing the dynamical 
processes occurring during the transition. The shift in 7„ is small enough for us to ignore 
as long as 7/;e <C 10^. Furthermore, as long as 7i;e <^ mp/me, the corresponding shift 
in the proton distribution can also be ignored. To summarize, it appears that both the 
negative and positive charges maintain their outward relativistic motion, but on balance 
the electrons are either shghtly faster or shghtly slower than the protons in order to account 
for the required current Uc to maintain Bip. As we shall see, these two situations result in 
significantly different 7-ray emissivities, and this difference may therefore constitute an 
important observational signature. 

In the absence of any resistive forces, the electric field would be quenched once 
the required current is produced by the time the particles reach zq. However, cooling 
processes (primarily from inverse Compton scatterings with the disk photons) above 
Zq continue to downshift the energy distribution of the electrons and protons (Melia & 
Konigl 1989). As such, this "Compton resistance" results in the generation of an electric 
field to maintain the current throughout the entire emission region. (We assume that 
the relativistic proton distribution does not change due to Compton upscatterings since 
by comparison with the electrons, the proton Compton scattering cross section, and 
hence the proton resistivity, is significantly smaller.) To determine the effect of these 
interactions on the electron distribution, we first calculate the response of a single electron 
at {z, r, 4>) moving with Lorentz factor 7 through the radiation field in the z direction. 
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Primed quantities denote values in the electron rest frame, whereas unprimed parameters 
pertain to the disk frame (hereafter referred to as the lab frame). The radiation field 
above the disk is characterized by the specific photon number density per solid angle 
nph{e,T[R]) — {2e'^/h^c^){e:K.p{e/kT[R]} — where e is the lab-frame photon energy 
and T[R] is the radius-dependent temperature. The electron moving through this field 
scatters dN photons to energies between Ss and Es + dSg and solid angles between ij,s4>s and 
[/Is + d/is] [<ps + d4>s] at a rate (per energy per solid angle) 



dN 



dt dSg djig d(f)s 



1 de 1 


'li R dR d^' 


1 ph J disk 





nph{e,T[R]) 



dUKN \ C(l - Pjj) 



df^sd<l>'sdesJ 7(1 -/3/^s) ' 



(5) 



where /5 = (1 - 7 ^) d = {z"^ + [r cos(0) - R cos($)]^ + [r sin(0) - R sin(<l>)] 



2N1/2 : 



is 



the distance between the electron and disk element R dR d^ located at {R, $), and ji — z/d 
is the cosine of the angle between the direction of the scattered photons originating from 
the disk at (i?, $) and the electron's direction of motion. (Note that is also the cosine of 
the angle between the direction of the photons incident on the electron and the direction 
normal to the disk.) The expression in square brackets represents the solid angle subtended 
by the disk element R dR d^ at (i?, $) with respect to the electron. The differential 
Klein-Nishina cross-section daKN/dn'gd(j)'gde'g is evaluated in the electron rest frame. Using 
the general expressions relating the lab and rest frame energies {e' — ej[l — P/i]) and angles 
{/i' — [ji — /3]/[l — P/j]', (p' — (f)), one finds the relation desd^sd(t>s/ de'gd^'gd(f)'g — 7(1 — (5^s)-i 
thereby allowing Equation (4) to be easily integrated over all scattered photon energies and 
solid angles to yield the single electron scattering rate. From this follows the 7-dependent 
mean-free path length for the electron: 



Imfpil] Z-i ^) — 



dN{j; z, r) 


-1 


f de f 


> RdRd^' 


c dt 




.J ph J disk 


d'' 



nph{e,T[R]) aKN {e') {I - (3pi) 



(6) 



As a result of a single scattering with a photon of energy e and direction /i, an electron 
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with initial Lorentz factor 7 attains a lower Lorentz factor 7^ in accordance with the energy 
conservation equation 7s = 7 — [^7^(1 — /9/x)]/[meC^ + £7(1 — P^)], where for computational 
ease, we have set /I'g and ji' to their average values (i.e., ii'g = jl' = 0). Thus, scattering 
that occurs within a distance dz reduces the number density of electrons with Lorentz 
factors between 7 and 7 + c?7 by 715(7; r)d'y dz/lmfpi'J', z, r). Since the particle number is 
conserved, this decrease corresponds to an increase in the number density of electrons with 
a Lorentz factor between 7^ and 7^ + djs of rieij; z, r)djs dz/lmfp{j', z, r). In the absence 
of a restoring force (i.e., if there were no electric field), the average value of the Lorentz 
factor (7) would undergo a gradual downward shift. As mentioned above, however, this 
downward shift in the electron density function must be balanced by the presence of an 
electric field which maintains the average Lorentz factor (7) of the distribution, in order to 
sustain the current associated with B^. Since the electric field shifts all electrons moving 
through the same distance by the same A7, we first calculate the change in the electron 
distribution function as a result of scattering along a path element dz and then linearly shift 
the degraded distribution function upward so as to maintain the value of (7). Repeating 
this procedure stepwise moving outward from zq provides a determination of ne(7; z, r) for 
all z. 

With the electron distribution thus specified, the 7-ray Compton emissivity can be 
determined by integrating Equation (4) over the entire scattering electron population. 
However, while a unidirectional motion was assumed above, the sensitivity of the photon 
production rate to the angle between the electron direction of motion and the scattered 
photon direction requires a more realistic treatment of the electron angular distribution. 
We leave the details of this development to future work. The fact that a strong correlation 
is observed between 7-ray emissivity and blazar activity suggests that 7-rays are emitted 
within a relatively narrow cone (see, e.g., von Montigny et al. 1995a). We here assume that 
this cone has an opening angle ~ 5°, and average the overall photon production rate over 
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1 > A*s > cos 5°. Thus, the rate at which photons are detected by an observer at a distance 
D (ignoring redshift effects) is given by the expression 

dtdSsdA (1 - cos5'')D2 j^^ J^.^ J ^^^o * ^ ' ' dtdSsd/j^sdcps ' 

(7) 

where azimuthal symmetry has been invoked to perform the integral over and average 
over 0s and the relationship between the scattered photon solid angle d^s d4>s and the 
detector area dA (i.e., dA — D'^d^sd(l)s) has been used. 

3. Results and Discussion 

(a) Particle Distributions and Characteristic Spectra 

In this section, we present the results of a calculation for a 10^ Mq black hole 
surrounded by a disk with T^ax — 5 x 10^ K and inner and outer radii Rin — 3Rg and 
Rout = 500Rg (where in this case Rg ^ 3 x 10^^ cm). Although the luminosity associated 
with this system (L fa 2 x 10^^ ergs s~^) corresponds to the weakest observed AGNs, 
qualitatively similar results are expected to follow from the modeling of more powerful 
sources, though this will require substantially more extensive computational resources. The 
active "shock" region is assumed to extend out to a radius Ra = 10 Rin, and the distance 
to the source is set to 100 Mpc. The efficiency rj is taken to be 10% and the transition zone 
boundary zq is placed at an arbitrary value of 5Rg. The dependence of our results on the 
actual value of zq is weaker than their dependence on the other unknowns in the problem, 
such as 7i;e, as long as zq <^ Rout- We have therefore chosen to concentrate our simulations 
on a parameter search corresponding to the dominant variables. We note, however, that the 
radiation field from the disk is anisotropic, and a careful treatment of the transition layer 
may produce some differences with the results shown here. 
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Before discussing these results, let us first analyze the photon spectrum produced by a 
power-law electron distribution with spectral index a injected at Zq that completely cools via 
inverse- Compton scattering with the disk photons (i.e., the well-studied scenario in which 
there is no re-energizing, e.g., due to an electric field). The resulting radiation spectrum is 
characterized by a power-law {dNph oc ej^) with index T — {a + 2)/2 for energies above the 
characteristic value Sch = ^kTmax (~ 10~^ MeV for the parameter space considered here). 
In practice, however, only electrons with 7 greater than a critical value jc (~ 100 for these 
conditions) cool appreciably (see, e.g., Dermer & Schlickeiser 1993; von Montigny et al. 
1995a; and references cited therein). As such, the distribution function of the outwardly 
moving electrons will evolve until it exhibits a break at ~ 7c. The corresponding photon 
spectrum, in turn, will exhibit a break at energy ~ ^c^ch (~ 10~^ MeV), below which the 
spectral index will be F = (a -|- l)/2. This effect produces a broken power-law spectrum 
with a shift in F between the X-rays and 7-rays of AF ^ 0.5 (see cooling spectra for 7; = 1 
in figures [3a] and [3b]; hereafter, 7; always refers to the electrons). 

This situation changes considerably in the presence of an electric field that acts above 
the transition zone to maintain the value of (7) for the electron distribution injected at Zq. 
In order to offset the energy lost by the high-energy electrons (which cool very efficiently 
via their interactions with the radiation field) , an electric field is established that shifts all 
of the electrons upward, including the low-energy part of the distribution. The electrons 
thus "bunch" up at (7), which in turn produces a bump in the photon spectrum at an 
energy (7)^£c/i- The overall spectral shape of the upscattered photons depends on how 
the value of (7) compares with 7,,. This point is illustrated in Figures 1 and 2. Figure la 
shows the evolution of the electron distribution function at ten different values of z for 
the model parameters a — 2 and 7^ = 1. The curves in figure la correspond to the values 
{z - Zo)/Rg = 0.0027,0.011,0.030,0.080,0.21,0.53,1.4,3.5,9.1, and 23, with the lowest 
value corresponding to the power-law curve and the subsequent values pertaining to the 
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curves evolving away from this initial distribution. The value of (7) for this evolution is 
13.8. As expected, the low-energy electrons are shifted upward and bunch near (7) while 
the high-energy electrons cool to 7c. The resulting photon spectrum is shown in Figure 2. 
Since (7) < 7c, the bump at (7)^ £ch (~ 10~^ MeV) is succeeded by a curve that softens 
to a power- law with spectral index F fa 2 for energies greater than 7c£c7i (0-1 MeV). The 
resulting shift in spectral index (i.e., ~ 0.5) between the X-ray and 7-ray parts of the 
spectrum is similar to that found in the absence of an electric field. 

In contrast to the previous case where 7; = 1, Figure lb shows the 
evolution of the electron distribution function for ten different values of z 
when a — 2 and 7^ = 100. The curves in figure lb correspond to the values 
{z - Zq) /Rg = 0.0033, 0.013, 0.033, 0.090, 0.23, 0.58, 1.5, 3.7, 9.3, and 23, with the lowest value 
corresponding to the power-law curve and the subsequent values pertaining to the curves 
evolving away from this initial distribution. In this case, (7) ^ 920. The greater value of (7) 
over 7c results in all electrons bunching near (7). The resulting photon spectrum is shown 
in Figure 2. Since (7) > 7c, the bump at (7)^ Ech (~ 10 MeV) is succeeded by a power-law 
with spectral index F 2. This spectrum is not well represented by a broken power- law. 
We note here that the energized distribution (7/ = 100) is injected with approximately 70 
times the energy of the 7; = 1 distribution and as such produces ^ 70 times more high 
energy photons. This point will be elaborated upon later in the discussion. 

(b) Comparison with the Data 

The high energy photon spectra (30 MeV to 30 GeV) of AGNs observed by EGRET 
are well-represented by a power- law with a spectral index ranging from 1.4 to 3.0. This 
spectrum is a natural consequence of the inverse Compton scattering of thermal disk 
photons by a cooling power-law distribution of relativistic electrons, though this would 
seem to require a range in the electron distribution index a between 0.8 and 4.0 which may 
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be broader than that produced by shock acceleration. As we have discussed above, the 
problem with this scenario arises from the fact that a purely cooling electron distribution 
with 7i = 1 produces a spectrum that is not consistent with many AGN X-ray observations 
in that it predicts a larger emission of X-rays than is observed (von Montigny et al. 1995a). 
In contrast, it is clear from Figures 1 and 2 that a low energy cut-off in the electron energy 
distribution function, together with an electric field that maintains the average Lorentz 
factor of the outflowing electrons in the presence of inverse- Compton scattering with the 
disk radiation, can significantly reduce the predicted UV and X-ray spectral components. 
To illustrate this point in the context of the observations, wc present in Figure 3a the 
flux density spectra for the model parameters o; = 2 and 7; = 1, 10, and 100, along with 
the multiwavelength data from PKS 0827-1-243 and in Figure 3b the flux density spectra 
for the model parameters a — 3 and 7^ = 1, 10, and 100, along with the multiwavelength 
data from PKS 1741-038 (taken from Figure 5 of von Montigny et al. 1995a). The dotted 
lines extending from the three different spectra on each figure are for cooling populations 
with corresponding values of 7^ at the transition zone. The different scales were chosen 
for easy comparison of the data (represented by the left hand scale) to the model results 
(represented by the right hand scale). We have not attempted to optimize the fit here, so 
the difference between the observed and calculated fluxes is presumably due to an imprecise 
choice of the black hole mass M, source distance D, and choice of 7;. 

It is evident from our results that the presence of an energizing region which shifts 
the electron population upward is sufficient to produce a spectrum devoid of X-rays. The 
subsequent acceleration does, however, produce an observable signature in the form of a 
bump around 10 MeV. As one can see from Figure 3, the main effect of this feature on the 
high-energy component of the spectrum is a turn-up toward lower 7-ray energies. There 
may in fact already be some evidence of this in several of the EGRET AGNs discussed by 
von Montigny et al. (1995a). 
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(c) Gamma-ray Faint Versus Gamma-ray Loud Sources 

A second important prediction of this model is that because the current associated with 
Brj) is due to an imbalance between the proton and electron outflow rates, the sign of the 
charge flux should be opposite on either side of the disk. So the electrons are re-energized in 
the transition zone on one side, but the Lorentz factor of their counterparts on the reverse 
side must be down-shifted. Since the apparent luminosity is sensitive to the value of 7^, 
AGNs observed from the energized side of the disk are clearly much more likely to have 
a high (or observable) 7-ray flux, compared with those observed from the backside. The 
expectation is that only about half of the superluminal radio sources should therefore be 
observable in 7-rays, with a possible additional correlation between the spectral break and 
distance. 

In their analysis of superluminal and strong flat-spectrum radio quasars that were 
not seen by EGRET, von Montigny et al. (1995b) showed that about 35 sources that 
would otherwise be expected to be strong 7-ray emitters were 7-quiet, compared to roughly 
40 — 50 that were detected (either strongly or marginally). The lack of evidence for any 
deflnitive characteristics distinguishing these two groups led them to suggest two primary 
reasons for the difference: 1) time variability, and 2) orientation and beam geometry effects. 
The scenario we have developed in this paper would point to a third reason, namely that 
the undetected sources are for the most part 7-quiet due to the uni-directional current flow 
associated with B^. As we have discussed above, the electrons are somewhat decelerated 
relative to the protons on one side of the disk, rendering them ineffective in Compton 
upscattering the ambient low-energy radiation. 

However, since the observed radio emission associated with the superluminal motion is 
produced many parsecs away from the central engine, we would not expect in this model to 
see a correlation between 7-faintness and a lack of longer-wavelength blazar activity. Since 
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most of the mechanical energy flux is carried by the protons, the energy budget available 
for broadband emission in the more extended outflow is roughly equal on both sides of the 
disk. So the i?^-induced current flow may offer an elegant explanation for the distinction 
between 7-loud and 7-faint blazars. 

Finally, on a related matter, one may argue using the curves in Figure 3 that the 
X-rays may also be produced by the same mechanism that accounts for the 7-radiation, 
if the X-ray portion of the spectrum rises as indicated by the theoretical curves in this 
waveband. Note that to fit both the 7-ray and X-ray data, we would need 7; pa 50 for PKS 
0827+243 and 7/ fa 100 for PKS 1741-038. Unfortunately the available data are not yet 
sufficient for us to distinguish between this situation and one in which the X-rays constitute 
a completely separate component. Future observations may identify the X-ray spectral 
index with the required precision for us resolve this issue. 
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4. Figure Captions 

Fig. la - Evolution of the electron distribution function for ten different values 
of z when a = 2 and 7^ = 1. The curves correspond to {z — ZQ)/Rg = 
0.0027,0.011,0.030,0.080,0.21,0.53,1.4,3.5,9.1, and 23, with the lowest value indicating 
the initial power-law and subsequent values pertaining to the curves evolving away from 
this. The units on the vertical scale are arbitrary. 

Fig. lb - Evolution of the electron distribution function for ten different values 
of z when a — 2 and 7; = 100. The curves correspond to {z — zo)/Rg — 
0.0033,0.013,0.033,0.090,0.23,0.58,1.5,3.7,9.3, and 23, with the lowest value indicating 
the initial power-law and subsequent values pertaining to the curves evolving away from 
this. 

Fig. 2 - The photon spectrum as seen by an observer at a distance of 100 Mpc for o; = 2 and 
the two indicated values of ji. These spectra are associated with the particle distributions 
depicted in Fig. 1. 

Fig. 3 - (a) Flux density spectrum as seen by an observer at a distance of 100 Mpc for 
a = 2 and the three indicated values of 7; (right-hand scale) . Multiwavelength observations 
of PKS 0827-1-243 are included for comparison (left hand scale) . The data are taken from 
von Montigny et al. (1995a). (b) The flux density spectrum as seen by an observer at 
a distance of 100 Mpc for a = 3 and the three indicated values of 7^ (right-hand scale). 
Multiwavelength observations of PKS 1741-038 are included for comparison (left hand 
scale). The data are taken from von Montigny et al. (1995a). 
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